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Abstract 

We consider a class of matrices with a specific structure that arises, 
among other examples, in dynamic models for biological regulation of 
enzyme synthesis |BJ. We first show that a stability condition given in 
[2j is in fact a necessary and sufficient condition for diagonal stability of 
this class of matrices. We then revisit a recent generalization of jS] to 
nonlinear systems given in , and recover the same stability condition 
using our diagonal stability result. Unlike the input-output based 
arguments employed in [3], our proof gives a procedure to construct a 
Lyapunov function. Finally we study static nonlinearities that appear 
in the feedback path, and give a stability condition that mimics the 
Popov criterion. 

Main Result 

The results of this note were triggered by the recent paper [S] and by several 
discussions with its author. We give our main diagonal stability result in 
Theorem ^ below, and present its implications for stability of a class of 
interconnected systems in the form of corollaries to this theorem. 
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Theorem 1 A matrix of the form 



A 



-1 ••• -7i 

72 -1 '■■ 

o 73 -l ••• ; 

; ■-. ■■. ■-. o 

o ••• o 7n -1 



7» > °) i 



(1) 



is diagonally stable; that is, it satisfies 

DA + A T D < 
for some diagonal matrix D > 0, if and only if 

7i . . .7 n < sec(7r/n) n . 



(2) 
□ 



The proof is given in an appendix. The results surveyed in [HI E] for 
diagonal stability of various classes of matrices do not encompass the specific 
structure exhibited by (JTJ) . In particular, the sign reversal for 71 in (JTJ) rules 
out the "M-matrix" condition, which is applicable when all off-diagonal terms 
are nonnegative. 

We now apply Theorem ^ to characterize the stability of the feedback 
interconnection in Figure^ When each block Hi is a first-order linear system 
with transfer function Hi(s) = 7i/(rjS + 1), ji > 0, Tj > 0, then a state- 
space representation of the interconnection would be obtained from the A- 
matrix in by multiplying its ith row by 1/rj for i = l,---,n. Since 
row multiplications by positive constants do not change diagonal stability, 
Theorem recovers the result in [H] which states that if (0) holds then the 
feedback interconnection in Figure ^ is asymptotically stable. Theorem ^ 
further shows that stability can be proven with a Lyapunov function V = 
x T Px in which P is diagonal. 

The linear result in [B] has been extended in [S] to the situation where 
H^s are not restricted to be linear, and instead, characterized by the output 
feedback passivity (OFP) jlj (a.k.a. output strict passivity [7j) property: 



- P < -\\yi\\ 2 + ii < Ui,yi > 



(3) 
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where || • || and < •, • > denote, respectively, the norm and inner product in the 
extended L 2 space, and (3 > represents the bias due to initial conditions. 
Using this property, [5] proves that the secant condition (j2J) insures stability 
of the feedback interconnection in Figure ^ 

Unlike the input-output proof given in [HJ, we now assume that a storage 
function Vi is available for each block in Figure ^ an d show that a weighted 
sum of these V^'s, 

V = £<W, (4) 

i=l 

where di > are chosen following the procedure below, is a Lyapunov func- 
tion for the closed-loop system. Indeed, a storage function verifying the OFP 
property (jSJ) satisfies 

Vi < -V* + liUiVi (5) 
which, when substituted in (jlj) along with the interconnection conditions 

ui = -y n , Ui = yi-x, i = 2,---n, 

results in 

V < -y T DAy (6) 

where A is as in , and D is a diagonal matrix comprising of the coefficients 
di in (£Q). It then follows from Theorem^ that positive g^'s that render the 
right-hand side of © negative definite indeed exist if (j2J) holds: 

Corollary 1 Consider the feedback interconnection in FigureU\ and let Ui, Xi 
and yi denote the input, state vector, and output of each block Hi. Suppose, 
further, there exist C 1 storage functions Vi(xi), satisfying ([3J) with ji > 
along the state trajectories of each block. Under these conditions, if HJ) holds 
then there exist di > 0, i = 1, • • • ,n, such that the Lyapunov function 
satisfies 

n 

V = Y,diV<-t\y\ 2 
i=i 

for some e > 0. □ 

Corollary ^ still holds when some of the blocks are static nonlinearities 
satisfying the sector condition 

< -y 2 i + JiUiyi, 7i > 0, (7) 
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Figure 1: Feedback interconnection for Corollary ^ 



rather than the dynamic property To see this we let X denote the subset 
of indices i which correspond to dynamic blocks Hi satisfying ©, and employ 
the Lyapunov function 

V = J2 d i v i- ( 8 ) 
For the static blocks, that is Hi, i ^ X, we note from (JZJ) that the sum 

n 

d i{-y 2 i + li U iVi) di>0 (9) 

i = 1 

is nonnegative and, hence, 

n n 

V < d iVi + J2 d i(-Vi + liUiVi) < J2 d i(~Vi + liUiVi) = -y T DAy 
iei % = l i=l 

% i X 

(10) 

as in Then, as in Corollary ^ condition (J2J) insures existence of a D > 
such that V < — e\y\ 2 for some e > 0. 



A Popov Criterion 

A special case of interest is the feedback interconnection in Figure 121 where 
H^ i = 1, ■ ■ ■ , n, are dynamic blocks as in (J5J), and the feedback nonlinearity 
ip(t,-) satisfies the sector property: 

0<Vn1>(t,y n )<K£, (11) 
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rewritten here as 

< -ip(t,y n ) 2 + KiP(t,y n )y n . (12) 

If we treat the feedback nonlinearity as a new block y n+ \ = if)(t, y n ), and note 
from (JT2J) that it satisfies (JJJ) with 7 n+1 = n, we obtain from Corollary [T] and 
the ensuing discussion the stability condition: 

K7i---7n < sec(7r/(n + l)) (n+1) . (13) 




Figure 2: The feedback interconnection for Corollary 



This condition, however, may be conservative because it does not exploit 
the static nature of the feedback nonlinearity. Indeed, using the Popov Cri- 
terion, the authors of 6\ obtained a relaxed condition in which n + 1 in 
the right-hand side of ()13|) is reduced to n when H^s are first-order linear 
blocks of the form Hi(s) = r y%/{riS + 1) and the feedback nonlinearity is 
time-invariant. 

To extend this result to the case where ifj's are OFP as in and not 
necessarily linear, we recall that the main premise of the Popov Criterion is 
that a time-invariant sector nonlinearity, when cascaded with a first-order, 
stable, linear block preserves its passivity properties. This means that, by 
only restricting H n to be linear, and combining it with the feedback non- 
linearity as in Figure 01 the relaxed sector condition of [H] holds even if 
Hi, ■ ■ ■ , H n _i are nonlinear: 

Corollary 2 Consider the feedback interconnection in Figure where Hi, 
i — 1, ■ ■ ■ , n — 1, satisfy (0) with C 1 storage functions Vi and 7« > ; H n is a 
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linear block with transfer function 

H n {s) = r n > 0, 7 „ > 0, (14) 

T n S + 1 

the feedback nonlinearity ib(-) is time-invariant and satisfies the sector prop- 
erty Under these assumptions, if 

«7i • • 'In < sec(7r/n) ri , (15) 

then there exists a Lyapunov function of the form 

V = J2 diVi + d n / ip(a)da, di>0, i = l,---,n, (16) 
i=i 1/0 

satisfying 

V < -e\(y ir --,y n - u ij(y n ))\ 2 
for some e > 0. □ 

Proof: Rather than treat H n and as separate blocks, we combine them 
as in Figure El 

Hn ■ \ Vn = ^%n), (17) 

and define 

V n = KT n \ tjj(a)da (18) 



which, from (fT7|) . satisfies 

V n = -Ky n ^(y n ) + ^ln^{yn)yn-l- (19) 

Because —Ky n ip(yn) < —''Piyn) 2 from (fT2|). we conclude 

V n < -1p{ynY + Kln1p{y n )yn~l = ~Vn + ^InVnVn-X-, (20) 

which shows that H n is OFP as in with % = ^ n n. The result then 
follows from Corollary ^ □ 

Corollary |21 can be further generalized to the situation where other non- 
linearities exist in between the blocks if,, i — 1, ■ ■ ■ , n, in Figure 121 If such 
a nonlinearity is in the sector [0, Kj+i], and is preceded by a linear block 
Hi(s) = t J^i ) then the two can be treated as a single block with % = K^iji, 
thus reducing n in the right-hand side of (J2J). 
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Figure 3: An equivalent representation of the feedback system in Figure El 



When H n is a linear block H n (s) 



In 



TnS+1 



, its series interconnection with the 



[0, k] sector nonlinearity ip(-) constitutes a dynamic block H n which satisfies 
© with 7„ = ivy n . 



The Shortage of Passivity in a Cascade of OFP Systems 

When the blocks Hi,...,H n each satisfy the OFP property (J5J), their cas- 
cade interconnection in Figure 0] inherits the sum of their phases and loses 
passivity. The following corollary to Theorem Q quantifies the "shortage" of 
passivity in such a cascade: 

Corollary 3 Consider the cascade interconnection in Figure^ If each block 
Hi satisfies (0) with a C 1 storage function Vi and ji > 0, then for any 

5 > 7i • • • 7n cos(7r/(n + l))^ 1 ), (21) 

the cascade admits a storage function of the form satisfying 

V < -e\y\ 2 + 5u 2 + uy n . (22) 

for some e > 0. □ 

Inequality is an input feedforward passivity (IFP) property jl] where 
the number 5 represents the gain with which a feedforward path, if added 
from u to y n in Figure El would achieve passivity. Corollary El thus shows 
that the cascade of OFP systems (JSJ in which 7i > represents an "excess" 
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Figure 4: The cascade interconnection for Corollary 01 



of passivity, satisfies the IFP property (|22j) with a "shortage" characterized 
by§IJ). 



Proof of Corollary 02 Using (j3J), ©, and substituting Ui = yi-i, i = 
2, ■ ■ ■ , n, we rewrite f!22D as 

n i 

di(~yf + liViu) + J2 dii-Vi + liViVi-i) + $(-u 2 - -uy n ) < -e\y\ 2 . (23) 

i=2 6 

To show that di > 0, i — 1, ■ • • , n, satisfying (j^j) indeed exist, we define 



.4 



; ■•. ■•. ■•. o 

o •■• 7n -1 
and note that the left-hand side of (|2"3|) is 



-1 ••• 

7i -1 "•■ 

72 -1 '•■ 



(24) 



[u y T ]DA 



u 

y 



(25) 



where D := diag {6, di, ■ ■ ■ , d n }. Because A is of the form (|TJ) with dimension 
(n + 1), an application of Theorem shows that a diagonal D rendering 
(j2SJ) negative definite exists if and only if (71 • • -7 n f) < sec(7r/(n + l))( n+1 ). 
Because this condition is satisfied when 5 is as in (|21j) , we conclude that such 
a D > exists and, thus, (j22J) holds. □ 
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APPENDIX: Proof of Theorem Q] 

Necessity follows because, as shown in [6J, if (J2J) fails then A is not Hurwitz. 
To prove that (J2J) is sufficient for diagonal stability, we define 



2 ' 



r := ( 7 i..-7n) 1/n >0 

A := diagjl, -2SL, ^ 3 
I r 

and note that 



- A^AA 

Thus, with the choice 
we get 



(26) 



,i+i72---7i 



i 

r 




i 



r 
••■ 



... o 

1 '•• 

r 



-l ' ' v 

^-l) n+1 r 



1 



„n— 1 



(27) 



D = A 



-2 



(28) 

DA + A T D = A- 1 (A~ 1 AA + AA T A- 1 )A- 1 (29) 

which means that DA + A T D < holds if the symmetric part of ([27)1. given 
by 

(30) 



-(-A _1 AA — AA T A _1 ), 

is positive definite. To show that this is indeed the case, we note that (j2] 
exhibits a circulant structure pQ when n is odd, and a skew-circulant structure 
when n is even. In particular, it admits the eigenvalue-eigenvector pairs 



A fe = 1 + re 1 n 



1 ri -i — k -i2^-k -i(n-l)22LfclT r. 

— lie « e « • • • e v ' n A; 
n 



n 



when n is odd; and 

? \ n. n. ' 



1 + re H ™ + »> 



1 

n ' 



n -*(-+— *) -«2(^+^fc) 



when n is even. Since, in either case, (J27J) is diagonalizable with the unitary 
matrix V = [v± ■ ■ ■ v n ] , the eigenvalues of the symmetric part (|3U|) coincide 
with the real parts of A^'s above. Finally, because 



min Re{l + re 1 « }= min i?e{l + re 

fc=l,---n fc=l,---n 



} 



r cos(7r/n), 



we conclude that if (J2J) holds, that is r < sec(7r/n), then all eigenvalues of 
(|50|) are positive and, hence, (|HU|) is positive definite and is negative 
definite. □ 
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